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Abstract 

When two solid surfaces are brought in contact, water vapor present in the ambient 
air may condense in the region of the contact to form a liquid bridge connecting 
the two surfaces : this is the so-called capillary condensation. This phenomenon 
has drastic consequences on the contact between solids, modifying the macroscopic 
adhesion and friction properties. In this paper, we present a survey of the work 
we have performed both experimentally and theoretically to understand the micro- 
scopic foundations of the kinetics of capillary condensation. From the theoretical 
point of view, we have computed the free energy barrier associated with the conden- 
sation of the liquid from the gas in a confined system. These calculations allow to 
understand the existence of very large hysteresis, which is often associated with cap- 
illary condensation. This results are compatible with experimental results obtained 
with a surface forces apparatus in a vapor atmosphere, showing a large hysteris of 
the surface energy of two parallel planes as a function of their distance. In the second 
part, we present some experiments on the influence of humidity on the avalanche 
angle of granular media. We show that the ageing in time of this avalanche angle 
can be explained by the slow kinetics of capillary condensation in a random confined 
geometry. 
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1 Introduction 



Molecules confined in narrow pores, with pore widths of a few molecular diam- 
eters, can exhibit a wide range of physical behavior. The introduction of wall 
forces, and the competition between fluid-wall and fluid-fluid forces, can lead 
to interesting surface driven phase changes, since for a small confinement the 
surface effects can be more important than the bulk effects [1]. Such effects can 
be observed in porous materials which have a large specific area. Porous ma- 
terials are involved in many physical, chemical or biological processes. Their 
adsorption properties are known to present a variety of behavior related to the 
texture of the porous matrix, which provides an experimental way to analyze 
the pore size distribution. Interpretation of adsorption isotherms in these ma- 
terials commonly involves a well known phenomenon, capillary condensation 
[1-3], which corresponds to the condensation of liquid bridges in the pores. 
More fundamentally, capillary condensation is a gas-liquid phase transition 
shifted by confinement. A basic model of confinement is provided by the slab 
geometry, for which the fluid is confined between two parallel planar solid 
walls. The classical theory of capillarity [3] predicts that in this geometry the 
liquid phase condenses when the substrate-liquid surface tension jsl is smaller 
than the substrate- vapor surface tension 7 5 y, and when the distance between 
the surfaces is lower than H c satisfying the Kelvin equation: 

Ap Ap ~ 2( 1S v ~ 1sl)/H c (1) 

Here, Ap = p^ — py is the difference between the bulk densities of the liquid 
and the gas phase, Ap, = p sat — p is the (positive) undersaturation and p sat is 
the chemical potential at bulk coexistence. If the vapor can be considered as 
an ideal gas, we have: Ap = kBT\n(p sat /p vap ), where ks is the Boltzmann's 
constant, T is the absolute temperature and p sa t/Pva P the saturated vapor 
pressure divided by the partial pressure of the vapor. Although the equilib- 
rium properties of this transition have motivated many experimental [4-7] and 
theoretical studies [8,3,9], capillary condensation presents remarkable dynam- 
ical features which are still to be explained. The most striking feature is the 
huge metastability of the coexisting phases, which contrasts with the bulk 
liquid-vapor transition. 

Since capillary condensation is a first order phase transition, one should be 
able to identify a critical nucleus and a corresponding free energy barrier away 
from the spinodal. For sufficiently small H, it can be shown that the liquid 
films coating the solid surfaces become unstable due to fluid-fluid interactions 
and grow to fill the slab. This has been carefully studied by several authors 
[7,11,12]). In this article we show that, as in the homogeneous nucleation case, 
the shape of the critical nucleus results from the balance between surface and 
volume contributions. The height of the activation barrier is so large that 
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it can induce a large metastability of the vapor phase. This first theoretical 
result is compared to experiments on capillary condensation in a surface forces 
apparatus. In the second part of this article, we will discuss the influence of 
this slow kinetics of capillary condensation on the mechanical properties of a 
granular material in a humid atmosphere. 



2 Homogeneous nucleation of a liquid phase between ideal surfaces 

2.1 Method 

Since capillary condensation occurs only in a confined geometry, the problem 
which arises in computing an energy barrier for the vapor/liquid transition is 
the validity of the macroscopic concepts of the classical theory of capillarity. 

To address this problem we use the following approach: 

i We use a Density Functional Theory (DFT) model for the fluid phase, tak- 
ing into account the long range interactions with the solid surfaces, and study 
the time evolution of a metastable confined vapor with a Langevin equation. 
We perform this study in a two dimensional geometry and determine the en- 
ergy barrier associated with the vapor/liquid transition as a function of the 
distance between the walls. 

ii- We compare these results with the prediction of the classical theory of cap- 
illarity in two dimensions. The classical theory gives the correct qualitative 
behavior and dependency of the energy barrier as a function of the confine- 
ment up to a numerical prefactor. 

iii- We then use the classical theory to calculate the energy barrier to conden- 
sation between two parallel plates in the three dimensional case, and discuss 
the order of magnitude obtained. 

2.2 Theoretical calculation of the activation energy in a 2D geometry 

We first use a mesoscopic Landau- Ginzburg model for the grand potential of 
the 2D system confined between two walls. In terms of the local density p(r), 
we write the "excess" part of the grand potential Vl ex = Vl + p sa tV, where p sat 
is the pressure of the system at coexistence, as: 

V ex = jdr^\Vp\ 2 + W(p) + (Ap + v ext (z))p} (2) 

In this equation, m is a phenomenological parameter allowing a simple treat- 
ment of inhomogeneous fluids; v ext (z) is the confining external potential, which 
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Fig. 1. Snapshots of the density p{r) of the fluid for different times 
t = 2, 6, 12, 18, 20, 28. If a point appears in white it means that the local density 
at this point is higher than 0.8. To obtain p(t), we make an average on the mean 
density in the pore on different realizations of the noise. 

we took for each wall as v ext (z) = —e(a/(Az + <r)) 3 , with Az the distance 
to the corresponding wall; e and o have the dimensions of an energy and a 
distance, W(p) can be interpreted as the negative of the "excess" pressure 
PsatP — f(p) — Psat, with f(p) the bulk free-energy density [13]. To allow a 
phase transition in this system, we assume a phenomenological double well 
form for W(p) : W(p) = a(p — pv) 2 (p — Pl) 2 , where a is a phenomenologi- 
cal parameter [14]. The system is then driven by a non-conserved Langevin 
equation for p: 

f^ + ,M (3) 

where T is a phenomenological friction coefficient and rj is a Gaussian noise 
field related to T through the fluctuation-dissipation relationship [16]. This 
time dependent Landau-Ginzburg model provides a good phenomenological 
description of the dynamics of the density field p(f, t) as soon as the dynamics 
of the nucleation is not limited by the crossing of the activation barrier. This 
model has been previously successfully applied to homogeneous nucleation 
[15]. The units of energy, length are such that a — e — 1. Time is in units 
of t = (Tea 2 ) -1 with F — |. In these units, we took m = 1.66, a = 3.33, 
pi = 1, pv = 0.1. Typical values of the chemical potential and temperature 
are A/i ~ 0.016, T ~ 0.06 (which is roughly half the critical temperature in 
this model). Periodic boundary conditions with periodicity L x were applied in 
the lateral direction. The simulated system is initially a gas state filling the 
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whole pore, and its evolution is described by equation (3). A typical evolution 
(see figure 1) of the mean density in the slit shows that: i) As expected [8], 
due to the long range nature of the external potential a thick liquid film of 
thickness i rapidly forms on both walls on a short time scale T\ (£ ~ 3.8a and 
T\ ~ 5to in our case) ii) Fluctuations of the interfaces around their mean value 
I induce after a while a sudden coalescence of the films. This second process 
has a characteristic time r. 
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Fig. 2. Logarithm of condensation time as a function of the "effective" width of 
the slab H — 3£ for fixed Afi = 0.016. The dashed line is the theoretical prediction 
ln(r) = ln(ro) + a(H — 3£) 3 / 2 . The two parameters ln(ro) and a have been obtained 
from a least-square fit of the data in a hi(r) versus (H — 3£) 3//2 plot. 

Studying the influence of the temperature on r, we have shown [17] that this 
time r obeys an Arrhenius law r = tq exp(Af2t //cgX"), where AQ^ is identified 
as the energy barrier for nucleation. In the 2D case, we have calculated [17] 
the energy barrier expected from the classical theory of capillarity: 

Afit = l(A /U A P7 Ly) 1/2 ^ 3/2 (4) 

The H dependence (A/i being fixed) on the activation time is plotted in fig. 2. 
As seen in fig. 2, the DFT model provides a good agreement with this classical 
prediction as far as the distances between the walls H is replaced by an effec- 
tive distance H — 3£ to take into account the presence of the wetting films. The 
prefactor a can be estimated from the data plotted in fig. 2, yielding a = 0.68, 
while the classical prediction gives a = 1.03 (here the liquid- vapor surface ten- 
sion -at finite temperature T = 0.06- has been computed from independent 
Monte-Carlo simulations of the model, yielding jlv = 0.8). The macroscopic 
theory thus gives a correct qualitative picture and a semi- quantitative agree- 
ment of the activation energy for capillary condensation. 
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Fig. 3. Picture of the critical nucleus for capillary condensation in three dimensions. 
R* represent the lateral extension of the critical nucleus (see text for details). The 
total curvature k of the meniscus is equal to ft = 1/R C = 2/H c . Note that in 3D, k 
is the sum of the in-plane and "axisymmetric" curvature. 

2.3 Nucleation between parallel plates in three dimensions 

We use then the classical theory of capillary to estimate the activation en- 
ergy for nucleating a liquid phase between two parallel plates. In the grand- 
canonical ensemble the critical nucleus corresponds to a saddle-point of the 
grand potential. We will consider the perfect wetting situation jsv — Isl + 
Jlv, although a generalization to the partial wetting case is straightforward. 
The grand potential of a pore partially filled with liquid may be written [8] 



where Vy (resp. Vi) is the volume of the gas (resp. liquid) phase and A$l, 
A S v and A LV respectively denote the total solid-liquid, solid-vapor and liquid- 
vapor surface area. The following expression is obtained for the "excess" grand 
potential, AQ t ot = — Fly, with Qy the grand-potential of the system filled 
with the gas phase only : 



where we have used pv — Pl — ApAfi. Within classical capillarity, long range 
fluid-fluid interactions are not taken into account, therefore the critical nu- 
cleus has to be a liquid bridge connecting the plates. One also expects this 
critical nucleus to exhibit rotational invariance, so that AQ t ot in eq. (6) is 
best parameterized in cylindrical coordinates (see fig. 3). In terms of p(z), the 
position of the LV interface, one obtains 



Q = -pyVy - p L V L + -fsvAsv + ISL^sl + 1lvA LV 



(5) 



tot — 



1lvA lv + ^lvAsl + A/iApV L 



(6) 
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where the index z denotes differentiation. Extremalization of the grand po- 
tential (7) leads to the usual condition of mechanical equilibrium, the Laplace 
equation, which relates the local curvature k to the pressure drop according to 
the Laplace law of capillarity: r )LV K — Ap ~ ApAp. This condition remains 
valid although the nucleus corresponds to a saddle point and not a minimum 
of the grand-potential. 

The main difference with bulk homogeneous nucleation comes from the pres- 
sure drop at the interface: here, the liquid pressure inside the meniscus is 
lower than the gas pressure since \i < jj, sat , so that the critical nucleus takes 
the form of a liquid bridge between the solid substrates instead of the spheri- 
cal shape encountered in bulk homogeneous nucleation. The previous Laplace 
equation is non-linear and cannot be solved analytically. From dimensional 
arguments however, one expects AQ tot = •jLvH^f(H/H c ), with f(x) a dimen- 
sionless function. The latter can be obtained from the numerical resolution 
of the Laplace equation, yielding the shape of the meniscus [18]. Numerical 
integration of eq. (7) then gives the corresponding free energy barrier. The re- 
sult for the energy barrier AQ^ is plotted in figure 4. As can be seen from the 
figure, a divergence of Af^ is obtained as the pore width H reaches H c . When 
the axisymetric extension of the bridge R* = p(-f-) is large compared to H, the 
negative (axisymmetric) contribution to the curvature is negligible and the L- 
V profile can be approximated by a semi-circular shape. This allows to obtain 
explicit expressions for the different contributions to AQ tot in eq. (6) as a func- 
tion of the extension R* of the bridge, namely V L = nR* 2 H - ^R*H 2 + f # 3 , 
A S l = 2ttR* 2 and A LV = n 2 R*H — -kH 2 . Maximization of AQ tot as a function 
of R* yields the following expression for the free energy barrier 

which does exhibit a divergence at H ~ H c — 2^ LV / ApAji. As shown in 
figure 4, this approximate expression is in very good agreement with the 
numerical calculation, even at small confinement H . Physically, an impor- 
tant consequence of the diverging energy barrier at H c is that the gas phase 
becomes extremely metastable: for water at 25°C, at a relative humidity of 
Pvap/Psat = 40%, we obtain H c ~ 2 nm and ^lvH 2 ~ 70k B T. This numeri- 
cal estimate shows that the energy barrier is always larger than the thermal 
energy of the system, except when H/H c <^ 1. 

2.4 Experimental evidence of the metastability 

Metastability effects in capillary condensation are very often observed in porous 
media, where they are responsible for hysteresis loops in adsorption/desorption 
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Fig. 4. Free energy barrier (in 3D) as a function of the normalized width of the pore, 
H/ H c . The solid line is computed by numerical integration of the Laplace equation. 
The points are obtained from the analytical expression, eq. (8). 

isotherms. In such hysteresis loops however it is not possible to tell which 
branch is the stable one. We have studied the capillary condensation of alkane 
vapors between smooth metal surfaces with a Surface Forces Apparatus [7]. 
In this type of experiments, the pressure of vapor is kept constant, and the 
distance between the surfaces (a sphere and a plane) is varied. The conden- 
sation of a liquid phase is revealed by a very strong adhesion force due to 
the capillary depression in the liquid phase. On the contrary, almost no in- 
teraction is measured when the surfaces are separated by the vapor. In these 
experiments, one observes the hysteresis loop of the force as a function of 
the distance between the surfaces associated with the metastability effects in 
the liquid/vapor transition. The particularity of SFA experiments is that it is 
possible to know which branch of the curve correspond to the lower energy 
state and therefore which phase is stable for a given confinement. This is due 
to the so-called Derjaguin approximation, which relates the force F measured 
between a sphere of large radius of curvature R, and a plane, to the interaction 
energy per unit area U of two parallel plates separated by the same distance 



The force curve shows then clearly that for a large range of values of the 
confinement H < H c the vapor phase is metastable. In the hysteresis range, 
the vapor never condenses over the time of the experiments-typically some 
hours-showing that the energy barrier to overcome in order to form a liquid 
bridge between the surfaces is actually quite large. 



H[2): 



F(H)/R = 4nU(H) 



(9) 
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Fig. 5. Force between a sphere of pyrex of radius R = 3.29 mm covered by a 
platinum layer and a plane of pyrex in the presence of a vapor of n-heptane. The 
arrows indicate the direction of the surfaces. H c is the critical distance at which 
there never is a meniscus and H s the distance at which a liquid bridge appears 
when the surfaces are brought into contact. 

3 Slow kinetics of capillary condensation: the aging of the friction 
coefficient 



3. 1 Humidity induced aging of the avalanche angle of a granular medium 



Since the condensation of liquid in a confined geometry can be hindered by 
high activation energy, one expects that capillary condensation processes may 
display slow kinetics. We discuss here in more detail the influence of humidity 
an the slow evolution in time of a macroscopic property: the avalanche angle 
of granular media. 

Solid friction properties of a large number of solid materials are well described 
by the Coulombs'law, which states that the minimum tangential force T that 
must be applied to a contact between solids in order to bring them into relative 
motion obeys the relation T = [x s N, where N is the normal load applied to 
the contact and /i s the static friction coefficient. Systematic studies on various 
material (for a review see Baumberger [19]) show that fi s increases slowly 
with the time during which the solid have been into contact at rest. This 
aging behavior of the static friction coefficient is often found to be logarithmic 
in time [20]:/i s = b + a log(t„,) 

More recently, the influence of humidity on this aging behavior has been 
brought into evidence. Dieterich [21] first showed in rock onto rock friction 
, that aging occurs only in humid atmosphere. In our group, we have studied 
more specifically the influence of humidity on aging of friction properties in 
granular media. For that purpose, we have measured the avalanche angle of 
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Fig. 6. Experimental setup. The granular material is placed in a cylindrical drum 
with an inner diameter and a length of 10 cm. The lateral faces of the cylinder 
are made of glass with an opening at the center to allow an exchange with the 
outer atmosphere. The drum is half filled with glass beads, whose diameter ranges 
between 125 and 160 microns. The relative humidity H = p vap /p S at is kept at a 
constant value by the method of the saturated aqueous solutions of inorganic salts 
[22]. 

a granular media consisting of micrometric glass beads in a controlled envi- 
ronment (see fig. 6). Using a rotating drum, the time t w during which the 
granular heap remains at rest can easily be varied. 



The evolutions of the maximum angle of stability as a function of the resting 
time t w , at different humidities are shown in figure 7. The data clearly show 
an aging behavior, since the angle of first avalanche systematically increases 
with the resting time. 

A first important point in this aging behavior is that it cannot be explained 
by an increase of the friction coefficient alone. Indeed, avalanche angle larger 
than 90° can be obtained for high humidities and large waiting time; in those 
cases the granular heap has enough internal cohesion to remain stuck to the 
upper part of internal wall of the drum . The effect of such a cohesive force 
can be included in Coulomb's analysis by adding an additional force F c to the 
normal component of the weight of a layer of glass beads at the surface of the 
heap [23]. This yields the following condition for an avalanche to occur: 

tanfl = /i s (l+ Fc - ) (10) 
\ mg cos v j 

where mg is the weight of the layer of glass beads undergoing the avalanche. 
Since the data show a linear dependence of the tangent of the avalanche an- 
gle with \og(t w )/ cos 8, we can conclude that there is a cohesion force in the 
granular media which increases as F c oc logt w . This dependency is observed 
over more than 4 decades of time (resting times range from 5 to 200 000 s). 
We do not observe any saturation of the angle at times as long as 200 000 s. 
The amplitude of this aging in characterized by the slope of the lines obtained 
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Fig. 7. The tangent of the maximum stability angle as a function of where t w is 
the resting time of the pile of glass beads. The lines are least-square fits of the 
experimental data, whose slope is identified with a. Different humidity are plotted 
on the same figure : H = 3%, H = 22% and H = 43%. 
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Fig. 8. Variation of the slope a(H) characterizing the aging behavior of the first 
avalanche angle with the relative humidity pvap/psat. 

is this representation. 

A very important parameter for this aging behavior is the humidity of the sur- 
rounding atmosphere. We do not observe any aging in dry atmosphere. Figure 
8 shows that this aging increases with the relative humidity: this increase is 
reversible. 



3.2 A simple model of thermally activated process for aging 



Although a lot can be learned from the perfectly flat slab geometry, the latter 
is certainly too idealized to account for the kinetics of adsorption in "real" 
experimental systems such as the previous ones. 9 
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Fig. 9. Schematic representation of a contact between rough surfaces. The fact that 
we consider a rough surface on a fiat surface does not restrict the generality of this 
approach. 

We show here that these logarithmic time dependence of the cohesion force 
between the surfaces may be understood by taking into account the influence 
of surface roughness on the dynamics of capillary condensation. Let us consider 
a simple model consisting of two surfaces facing each other and rough at the 
nanometric scale, as depicted on fig. 9. As emphasized in the introduction 
capillary condensation typically occurs in pores of nanometric size. We thus 
have to consider the roughness of the surfaces at the nanometer level. Here 
again we shall stay at a macroscopic description, and focus on a qualitative 
picture of the influence of roughness on the transition mechanism. Without 
loss of generality, one may consider that one of the walls is perfectly flat. When 
roughness is present, there is a broad range of gaps between the surfaces. In 
particular, there are regions where the two surfaces are in close contact. In 
such regions, condensation should take place on a very short time-scale. Thus 
at "early times" , one has to consider a set of wetted islands, which we shall 
consider as independent one of the other. Once these islands have formed, they 
should grow up to a point where the distance between the surfaces is equal 
to H c , so that a meniscus of radius R c = H c /2cos8 forms at the liquid- vapor 
interface, allowing for mechanical equilibrium. 

In doing so however, the wetted area has to overcome unfavorable regions 
where the distance between the two surfaces is larger then H c . Let us consider 
a specific jump over such a "defect" , as idealized in fig. 10. We denote by hd 
the "averaged" gap inside the defect (hd > H c ), and by ad its area. The free 
energy cost for the liquid bridge to overcome this defect is approximatively 
given by 



where Vd is the excess volume of the defect, va = a d (e d — H c ). We can thus 
estimate the time to overcome the defect as 



One may expect the defects to have a broad distribution of excess volume Vd, 
so that the activation times r are very widely distributed. After a time t w , only 



Aft 1 " ~ a d (A/iAp e d - 2(<y S v ~ Isl)) 
= v d A/iAp 





(12) 
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the defects with an activation time r smaller than t have been filled by the 
liquid phase. Using eq. 11 and 12 these have an excess volume Vd which verifies 
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Fig. 10. Schematic representation of a defect. Vd is the volume of the region where 
the height hd > H c . 

Vd < Vdmaxitw) = k B T (ApAp)' 1 \n(t w / r ) . The number of filled defects at a 
time t is then typically N(t) = Vd ma x(t w )/vo where vq is the typical width of 
the distribution of excess volume of the defects. Now, once a liquid bridge 
has bypassed a defect, the area wetted by the liquid increases by some typical 
(roughness dependent) a^. The time wetted area grows in time as: 



A{t w ) ~ N{t w )a d 



ln(^) (13) 



Ap/(k B T)Apv \t 

Assuming a typical radius of curvature A of the asperities, we find if A ^> r K , 
that ad ~ Ar^. If we now assume that the vapor is an ideal gas, then A/x = 
k B T\ia(p sat /p vap ). This give a cohesion force between two beads: 



F c (t w ) = ^A(t w ) 
r K 

" LlA ln(t w /r ) (14) 



Apv \n(l/H) 

This model gives the good dependency of the cohesion force with the resting 
time t w in qualitative agreement with the experiments. It is not possible to 
test directly this prediction since vq and A are not easily measured but we can 
check the dependency of a with the relative humidity. However we can check 
the internal coherence of this mechanism by estimating the order of magnitude 
of the cohesion force needed to induce aging of the avalanche angle. Using 
equation 10, the increase of the cohesion force during the time t w must be: 

. „ mq cosO „, . 

AF C = — Atan0(O (15) 

The numerical value of DeltaF c is thus of the order of the weight of a bead. 
On the other hand, the maximum value of the capillary force between two 
beads in contact is obtained when the beads are ideal sphere (no roughness): 
F c = 27i'y L vR = 3.4 x 10~ 5 N for beads of radius R = 75 pm in water 
vapor humidity. This is four order of magnitude larger than the weight of a 
bead:4.4 x 10" 9 N. 
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Therefore one see that a slow kinetics of capillary condensation between rough 
surfaces, leads to a logarithmic growth in time of the capillary force and can 
actually be responsible for the aging of the avalanche angle of granular me- 
dia, as well as more generally for the humidity induce aging of the friction 
properties of surfaces. 
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